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^^ , The structure of states of the perturbed p-spin spherical spin-glass is analyzed. At low enough 

^— ^ ■ free energy metastable states have a supersymmetric structure, while at higher free energies the 

^^ ' supersymmetry is broken. The transition between the supersymmetric and the supersymmetry- 

(~| ' breaking phase is triggered by a change in the stability of states. 
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^_^ . It is well known that mean- field spin-glasses can be divided into two broad categories. On the one hand there are sys- 

P5 ' terns with a low temperature thermodynamic phase described by full (or continuous) breaking of the replica symmetry 

^ , (FRSB). In these systems the static transition temperature Tg coincides with the temperature Td where the relaxation 

^ ' time diverges. Moreover, the dynamical correlation function exhibits single-step relaxation, either exponential (high 

'I"! , temperature), or power-law (low temperature). The paradigm of these spin-glasses is the Sherrington-Kirkpatrick 

. ■ (SK) model [ij. On the other hand, there are systems where the static and dynamic transitions occur at different 

jrt [ temperatures, namely Tg < Td- In these systems thermodynamics is described by a one step breaking of the replica 

r^ . symmetry (IRSB), and the dynamical correlation function develops, at low temperatures, a two-steps relaxation, 

i' ' typical of structural glasses. The main representative of this second class is the p-spin spherical model f3|. In mean- 

'■^ . field spin-glasses, it seems that a given type of static behaviour is always accompanied by a given type of dynamical 

^ ' behaviour. No mixed cases are known. 

O ' In disordered systems statics and dynamics are linked to each other by the structure of states: lowest lying states, 

i.e. ground states, dominate the statics, whereas the dynamics is ruled by metastable states with higher free energy. 
It is therefore reasonable that the origin of the two different behaviours described above, and thus of the connection 
between statics and dynamics, has to be sought in the structures of states of different systems. It is common wisdom, 
^ ' for example, that metastable states play a much greater role in IRSB than FRSB spin-glasses, and that for this reason 
dynamics is so much different from thermodynamics in IRSB systems. However, it would be important to have a 
mathematical formalization of the different structure of states in the two classes. Recent studies 0; S B Ifl S H S 
^ ! EIHIIl suggest that the Becchi-Rouet-Stora-Tyutin (BRST) supersymmetry jH, El IH 111 113 may provide such 
f^ ■ a formalization. 

^^ [ The BRST supersymmetry is a transformation which leaves invariant the action which one has to define in order 

^D to study the structure of metastable states in spin-glasses 0j|j;[l3) and in particular when the complexity, i.e. the 

entropy of the states, is computed. As always in physics, the symmetry of the action may or may not be shared 

by the actual physical states. More specifically, in spin-glasses the complexity is equal to the action evaluated in 

its stationary point: there are systems where the physical (stable) stationary point is symmetric, and other systems 

where it is not. In this last case the supersymmetry is spontaneously broken. After one year of debate, it seems now 

clear that the structure of states of the SK model breaks the supersymmetry [l^, uM i while in the p-spin spherical 

Q ' model the supersymmetry gives the correct set of physical states |3, ■ This fact is in agreement with the other known 

O ■ physical features of these systems. Supersymmetry breaking occurs because metastable states are marginal and very 

unstable under external perturbations [iflll3jll3> whereas well shaped states, with finite second derivative and large 

barriers surrounding them, preserve the supersymmetry. This is exactly what we expect from the different dynamical 

behaviours and role of metastable states, in the SK and p-spin models. 

Classifying spin-glasses according to whether or not the supersymmetry is spontaneously broken seems therefore a 
sound criterion to distinguish different structures of states, and thus different static and dynamic behaviours. In order 
to understand the physical differences of these two cases in a deeper way, it would be desirable to find a model where 
both structures of states are present, and to study its analogies with the two archetypes described above, the SK 
and the p-spin model. In this paper we introduce such a model and show that supersymmetric and supersymmetry- 
breaking states can coexist in the same system, albeit at different free energies. We will argue that such a situation 
may be rather general in spin-glasses, and that the SK and p-spin models correspond to the two extreme cases: the 
transition from supersymmetric to supersymmetry-breaking states takes place at the ground state in the SK model, 
whereas the transition occurs at the highest free energy in the p-spin spherical model. 

Before starting our analysis, let us make a general remark. There is a dynamical form of the supersymmetry, 
which arises when studying the Langevin dynamics of statistical models [17. 19| . The important difference with its 
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static counterpart described here, is that the dynamical supersymmetry is always broken at low temperatures, in 
all models. In other words, it is an intrinsic consequence of ergodicity breaking, linked to the break-down of the 
dynamical fluctuation-dissipation theorem [l^. On the other hand, the static supersymmetry we are discussing may 
or may not be broken, according to the particular structure of states displayed by each given model. 
The Hamiltonian of the model we want to study is given by, 

JV N 

H=- ^ Jt,...ip (J^^■■■(J^p- £ Yl Kl^-h ^l^ . . . Gl^ , (1) 

ii<---<jp Ii<-<It 

where the spins Ui are real variables subject to the spherical constraint ^^ tr^ — N, and the Gaussian random 
couplings Jii...i and Ki-^,,,i^ have variance p\/2NP~^ and r\/2N^~^, respectively. We are therefore perturbing the 
p-spin spherical model with an extra r-body interaction term. It is known that p + r spherical models may display 
a nontrivial thermodynamic behaviour when p > 3 and r = 2: in that case there is a transition between a IRSB 
thermodynamic phase (low e), to a FRSB phase (large e) [23. On the contrary, if both p and r are strictly larger 
than two, we expect the model to have a normal IRSB thermodynamic behaviour. This is the case we will analyze. 
What we will show is that, even in this apparently simple case, as soon as e 7^ the structure of metastable states 
becomes much richer than the unperturbed p-spin model. In the following, all plots are done for p — 3 and r = 4. 

Metastable states are the local minima of the mean-field free energy, that is the Thouless- Anderson-Palmer (TAP) 
free energy |2j. The first thing to do is therefore to work out the TAP free energy for this model. This can be done 
by following the method used by Rieger in |23|. We will work with e <C 1, and thus keep only the lowest order terms 
in e. The TAP free energy density is 23], 

B ^ ^ 1 

PfTAP^ ^ ]V ^ -'^ii-.-jp »7i,, ...m,p - e— ^ Ki^,„i^^mi^...mi^~ - log(l - g) 

ii<---<ip h<-<lr 

- ^[{P- 1)<Z^ - pq"-' + 1] - e' Y [(^ " ^^^^ " '''^'"' + ^] ' (2) 

where m^ = (cr^) are the local magnetizations, and q is the self-overlap of a state, q = J2i 'mf/N . We had to keep 
the O(e^) term since, after averaging over Ki-^,,,i^, also the linear term will become O(e^). To compute the TAP 
complexity we closely follow the standard method, as can be found, for example, in |j|. The number of TAP states 
with free energy density / is. 



mf)^ Jl[dm, 6[d,PfTApini)] det[d,dj (3 frApim)] S[(3fTAp{m) 



' PI] , (3) 



where, as usual, the modulus of the determinant has been dropped (for a discussion of this point see [laliSliSl ^^^ 
1^). We give an exponential representation for the two delta functions introducing the Lagrange multipliers Xi and 
u, and of the determinant, with the anti-commuting Grassmann vectors {"0^, tpi} H- ^^ ^^^^ ^^Y ^"^ '^^^ write. 



JV{f)= I VmVxVipVil^du e^("'^>''^'"'/) , (4) 

where the action S is given by, 

S{m,x,tp,-ip,u; f) ^ p'^XidifrApim) + (3'Y^ril)jdidjjTAp{m) + PufTAp{rn) - Puf . (5) 
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We work at the annealed level, which is exact for e = 0. After averaging the number M{f) over the disorder J^j.. 
and Ki-^,,,1^, the integrals in 'DmVx'D^p'D'ijj can be performed exactly, provided that we introduce the parameters. 



T= —^tp^ilj, B = —^x^m, q= —Ym.rrii . (6) 

i i i 

After some algebra we obtain, 

Wj) = I dTdB dq du e^^(^^S'9'"! /) , (7) 




FIG. 1: Unperturbed p-spin spherical model: the complexity as a function of the free energy density. The stability parameter 
Xp is zero at the threshold / = fth- The temperature is equal to the dynamical temperature Td = 0.612 of the unperturbed 
model. 



where the effective action S given by |23| . 

S = f3u[g{q) + e^h{q)-f] + {B^-T 
l.o _, 1 
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and where we used the following definitions 
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From (UJ we have that in the thermodynamic Umit N —^ oo the complexity is given by the effective action evaluated 
in the saddle point values of the variational parameters T, B, q, u, 



j:{f)^-\ogAf{f)^ S{T,B,q,u; f) 



saddle point 



(13) 



The saddle point equations are not too difficult to solve. In fact, the parameters T and B can be worked out explicitly 
as functions of q. On the other hand, the equation for q itself. 



dS{q,u;f) 



0, 



(14) 



is highly nonlinear and must be solved numerically. The complexity is naturally a function of the free energy density 
/, so in principle we have to minimize S also with respect to u, such that / remains the only free variable. However, 
/ and u are conjugated variables (this is clear from the form of the action in |SJ|). This means that we can either let 
/ free, work out u = u{f) and study S(/), or let u be the independent variable, such that / = f{u) and study I](u). 
It turns out that in general E is a smoother function of u than /, so we will mostly use this last representation. 

Action |(SJ) is invariant under the generalized Becchi-Rouet-Stora-Tyutin (BRST) supcrsymmetry [isL ll^-|a. 1^. 
This symmetry generates the following two Ward identities |3, |j| , 



{m ■ x) — —(ip ■ tp) 



{x ■ x) = {u ip ■ 4') ■ 



(15) 
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FIG. 2: Unperturbed p-spin spherical model. Complexity as a function of u. 
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FIG. 3: Unperturbed p-spin spherical model. Upper panel: the stability parameter Xp. Lower panel: the supersymmetry 
parameter ^. 



The first Ward identity immediately translates into an equation for the parameters T and B. From equations © we 
have, 



BSRT => B + T = . 



(16) 



Given that both T and B are functions of q, imposing this relation is equivalent to impose an extra equation for 
q, which has to be satisfied together with (|14|l . By doing this we select the supersymmetric solution of the saddle 
point equations. It is straightforward to prove that the second Ward identity in (|15|l gives a relation for q which is 
equivalent to (|16|l . 

The saddle point equations for the complexity may, in general, have many solutions. Some of these solutions 
are supersymmetric (SS), and therefore satisfy the Ward identities H15|l . whereas others are supersymmetry-breaking 
(SSB) solutions. These last solutions violate the Ward identities. Physically, different solutions correspond to different 
structures of metastable states, so selecting the correct solution, whether supersymmetric or not, is a key point. As 
we have seen, from the studies performed so far, it seems that there are two classes of models. In IRSB systems 
(as the unperturbed p-spin spherical model) the supersymmetry is unbroken 0; 01 j so one can directly impose the 
Ward identities (|15|l to simplify the saddle point equations and find the physical complexity. In the SK model (and 
presumably in other FRSB systems), the supersymmetric saddle point is unstable, and physical metastable states are 
described by the SSB saddle point. In this case imposing ifTB)) gives the wrong result. In all these cases, the main way 
to discriminate the physical from the unphysical solution is to test the stability of the states described by the various 
solutions. In particular, a very important quantity is the parameter introduced by Plefka in p4|. 
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In order to have stable states, with a physical susceptibility, we must have, 



Xp >0 



(18) 
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FIG. 4: Perturbed p-spin spherical model (e = 0.45). The complexity as a function of u. Full line: the SS solution. Dashed 
line: the SSB solution. Open circles indicate the physical solution. T = 0.6Td. 
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FIG. 5: Perturbed p-spin spherical model (e — 0.45). Upper panel: the stability parameter Xp. Lower panel: the supersymmetry 
parameter (,. T — 0.6Td. 



We shall use the stability condition H18|l to classify physical and unphysical solutions of the saddle point equation 
of our model. As a starting point we briefly discuss the case e = 0, i.e. the unperturbed p-spin spherical model |25| . 
As we have said, in this case the physical complexity is the one corresponding to the supersymmetric solution of the 
saddle point equations. Such solution can easily be found by imposing (|16|) . The complexity E as a function of the 
free energy density / (Fig.l) is zero at the ground state /o, and it is defined up to a threshold value fth, beyond which 
the saddle-point equations have no longer solution. Threshold states are marginal, i.e. they have a zero mode in the 
second derivative of the free energy, and Xp = (whereas Xp > for all states with / < fth) ■ 

The fact that threshold states are marginal suggests, in analogy with the SK model [lOllll]|. that they may have 
a supersymmetry-breaking nature. If we study the complexity as a function of the parameter u (Fig. 2), we find that 
I](m) increases with m up to a value Uth — u(fth), and beyond that point it remains constant, i.e. S(m) = "Eth for each 
u > Uth- Note that above Uth also the free energy is constant, f{u) = fth for u > uth, and thus the fiat branch of S(w) 
is concentrated in the threshold point of S(/) in Fig. 1. The flat branch of the curve has Xp = (Fig. 3), meaning, 
as we have already said, that threshold states are only marginally stable. To assess the supersymmetric properties of 
this curve, we can introduce the following parameter, 



i = T + B 



(19) 



From (|16|l we have that ^ = on a SS saddle point, and ^ 7^ on a SSB saddle point. We see from Fig. 3 that the fiat 
branch of S(u) has ^ 7^ and thus breaks the supersymmetry, while the growing part of the curve is SS. It should be 
noted that the fiat branch continues up to u = 0, which corresponds to putting no constraint on the free energy, i.e. 
to counting the total number of states. In other words, in the p-spin spherical model, the total, that is the maximum, 
complexity Y^th can be found either from the SS solution at u = uth, or from the SSB solution at u = 0. The situation 
in the unperturbed case seems therefore quite strange: if we use / as a natural variable, we have states up to / = fth, 
and these are all supersymmetric and stable. Only threshold states are marginally stable. On the other hand, using 
M as a variable, we find a degenerate interval between uth and 0, where the complexity and the free energy remain 



stuck at their threshold levels, fth and T,th, and where the supersymmetry is broken. This situation becomes clearer 
once we perturb the system. 

When we switch on the perturbation e we recognize that the e = case is just the peculiar (and rather pathological) 
limit of a more general physical situation. We first analyze S as a function of u (Fig. 4). There are two solutions 
of the saddle point equations for each value of u: the SS solution, and the SSB one. These two branches touch at 
u — uth{^)- Below this point, the parameter Xp is negative for the SSB solution, and positive for the SS one. On the 
other hand, Xp changes sign for both solutions at uth^ and thus the physical stability of the two solutions switches at 
Uth- As a result SS states are the physical one below uth, whereas SSB states become physical above uth- As it can 
be seen from the plot, the SSB branch is very weakly dependent on u when e is small, and in the limit e = it is just 
flat, as in Fig. 2. 

Therefore, threshold states are, in the general case, at the same time the highest SS and the lowest SSB states. 
This is in agreement with the marginal nature of these states, as it was already well clear in the unperturbed case. 
It is interesting to note that the SSB solution in its stable phase u > Uth is completely analogous to the SK model, 
suggesting that the same connection holds also for the physical properties of the two models in this phase. The fact 
that the SK model has a FRSB thermodynamics, while the present one is a strictly IRSB model, is simply due to the 
fact that for the present model, 

nuth) > , (20) 

such that the SSB states do not contribute to the statics. The statics of the SK model is FRSB simply because the 
threshold states in that model coincide with the ground states, i.e. 

E(ut,0 = SK . (21) 

On the other hand, in the unperturbed p-spin model we have the opposite extreme behaviour: there are no states 
whatsoever above the threshold, and therefore no SK phase, not even in the metastable spectrum of states. 

The crossover between a SS and a SSB phase can also be analyzed by looking at the behaviour of E as a function 
of /. Let us consider first the SS solution: the complexity grows from /o where it is zero, up to a value fth = f{uth), 
where u = uth- Beyond uth both E and / decreases, giving rise to a new branch of the SS curve E(/). This backward 
branch is SS, but unphysical, since Xp < 0. The SSB solution takes over exactly at fth, providing a continuation of 
E(/) beyond the threshold. This state of affairs is reproduced in Fig. 6, where we plot E in a neighborhood of fth- 
Therefore, unlike in the unperturbed p-spin spherical model, we have some physical states also above the threshold, 
and these states have a SSB structure. Note that in this context fth is not the largest free energy density of metastable 
states (as normally in the literature), but it is the largest free energy density of supersymmetric metastable states. 

Our result suggests that the most general structure of metastable states in spin-glasses is the following. At low 
free energies lie SS states: these are non-marginal minima of the free energy, well separated from each other by large 
barriers. The stability of these states decreases with increasing energy: the softest mode of their second derivative 
matrix is smaller the larger /. At the threshold energy SS states become unstable, and above this point SSB states 
are the physical ones. Presumably, as in the SK model, SSB states are marginal, in the sense that they have a zero 
mode leading out of the states, and therefore are not well separated by large free energy barriers. 

This crossover in free energy between SS and SSB states has important consequences on both the statical and 
dynamical behaviour of a system. As we have said above, when the transition free energy, i.e. the threshold fth 
is larger than the ground state free energy /o, SS states dominate the thermodynamic properties, and the model is 
statically IRSB. If, on the other hand, fth is equal to /o, the model has a FRSB statics. Note that the difference, 

A = fth-fo, (22) 

may depend on the external parameters, as the temperature, the field or any other. In the present model, by changing 
the temperature, A changes, but it always remains strictly positive, so that the model is always statically IRSB. 
However, for other systems the situation may be different. In particular in the Ising p-spin model, there is a transition 
between a IRSB and a FRSB phase at the Gardner temperature Tq 0. Moreover, it has been proved in [23| that in 
this model a crossover occurs from IRSB to FRSB metastable states in the complexity. Thus, what we expect is that 
in such a model the supersymmetric threshold becomes equal to the ground state at Tq, i.e. A(Tg) = 0. The study 
of the supersymmetric properties of the Ising p-spin model, will appear in [23 . Preliminary results are in agreement 
with what we have shown here. 

From a dynamical point of view, it is reasonable to assume that ergodicity is broken by the highest non-marginal, 
i.e. supersymmetric, states. Therefore, we expect that the dynamical transition is triggered by threshold states. In 
systems where fth > fo, as the spherical p-spin model, we will have Td > Ts, whereas in systems where fth ~ fa, as in 
the SK model, Td = Tg. We would hkc to remark that the high free energy transition between SS and SSB metastable 
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FIG. 6: Perturbed p-spin spherical model (e — 0.45). Complexity as a function of the free energy density, around / = fth- 
Inset: the fuU curve E(/). T = O.QTa 

states, and its consequences on the static and dynamical behaviour of spin-glasses, was for the first time discussed 
and understood in 27], albeit not in a supersymmetric context. In that study, the same phenomenon was interpreted 
as a IRSB to FRSB transition in the spectrum of states. The physics, however, is the same. 

In this paper we have analyzed a spin-glass model where a crossover takes place between low energies supersymmetric 
states, and higher energies supersymmetry-breaking states. The model is always solved, at the static level, by a IRSB 
thermodynamic ansatz, although SSB states have the same physical structure as in FRSB systems. This is due to 
the fact that the crossover always occurs at energies higher than the ground state. Our study and the results of [23 
suggest that this may be the generic scenario, and that different static and dynamic behaviours may just depend on 
whether the transition free energy, i.e. the threshold, is larger or equal than the ground state free energy. 
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